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ABSTRACT 



Short and long distance contributions to the exclusive B-decays into various 
K- resonances and dileptons, i.e. B — > K l ti{t = e,fj,,v), are examined. The 
heavy quark effective theory has been used to calculate the hadronic matrix 
elements. Substantial branching fractions are obtained for the dileptonic 
B-decays into some higher excited states of K-mesons. The long distance 
(resonance) contributions to these exclusive rare B-decay modes dominate 
the short distance contributions mostly by two orders of magnitude. It is 
pointed out that, excluding the resonance contributions, the P-wave chan- 
nels are dominant, accounting for about 50% of the inclusive B — > X s 
branching fraction. 



1 Introduction 



Rare B-decays have been in the focus of a lot of theoretical attentions . This is due 
to the extensive amount of information on Standard Model (SM) that can be extracted 
from these processes. The rare decays proceed through flavor changing neutral current 
(FCNC) vertices that are absent at the tree level, therefore, they are a good probe of the 
SM at the quantum (loop) level. On the other hand, radiative B-decays are sensitive to 
quark mixing angles Vtd, Vt. s and Vf&, hence their measurements yield valuable information 
on CKM matrix elements, and consequently, shed some light on the CP violation in SM. 
Also, reducing the uncertainties of the theoretical calculation of these processes within 
the SM, and comparing them with more precise measurements available at future B- 
factories, may be the most promising probe of the new physics in the short term. 

In this paper, we focus on the rare B-decays, B — > K % tt{t = e,/i, v\ where K % are 
various resonances of K- meson. We calculate both, the short-distance (SD), and the 
long-distance (LD) contributions to these processes. The motivation for this work is the 
following: Due to the large phase space available to rare B-decays, exclusive processes 
with excited K-mesons in the final state might have substantial branching fractions. In 
fact, in reference 0, it is shown that for the rare B-decays B — > K l j, which proceed 
through SD penguin operator, some higher resonances of K-meson have 2-3 times larger 
branching ratios than the ground state. On the other hand, in reference ||, we indicated 
that the LD contributions to these processes are very significant, approximately 20% of 
the SD amplitude. Consequently, we feel that the same investigation should be carried 
out on dileptonic rare B-decays. 

A main source of uncertainty, in all these estimates, has always been the evalu- 
ation of hadronic matrix elements (HME) for specific exclusive decays. This would 
involve the long range nonperturbative QCD effects which renders the calculations 
model-dependence. We use the heavy quark effective theory (HQET)||, as a convenient 
ansatz in formulating the HME. This results in an enormous simplification, especially 
when decay rates to higher excited K-meson states are estimated. With a plausible 
assumption that b quark is static within the B meson, one can show that these HME 
are the same as those appearing in B — ► K l tfj(ifj') decays. Therefore, one can adjust 
the universal form factor, such that a HQET formulation of the HME fits the exper- 
imental data available for B — ► T^W^^H]. j n f ac ^ our predicted value for the ratio 
R = T(B -> K*^)/T(B -> X,j), obtained this way (see eqn. (22) in Ref. 0), is in 
good agreement with the recent CLEO results |J. In a follow up paper, we showed, by 
adopting a more general model for the universal form factors, one can also estimate the 
non-leptonic decays to higher K-meson resonances ||. These predictions are yet to be 
tested experimentally. 
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2 Differential decay rates 



We start with the effective Lagrangian for b — > , which includes both, short and 

long distance contributions f7|, § 

L e// = ^ f^r) V^Vt^AsL^blLn + BsLfMRH + 2m b s 2 w FsT fl b£^£), (1) 

where 

L M = 7 M (1 -75), ^ = 7^(1 + 75), 

and 

Tn = -icr flu (l + 'j 5 )q u /q 2 . 

Vij are the Cabibbo-Kobayashi-Maskawa matrix elements, = sin 2 9 w ~ 0.23 (6^ is 
the weak angle), Gf is the Fermi constant and q is the total momentum of the final 
pair. 

This effective Lagrangian is obtained by integrating out the heavy degrees of freedom, 
ie. the top quark, W and Z bosons, at the scale /1 = My/. Using the renormalization 
group equations to scale down to a subtraction point comparable to the light masses, 
ie. /1 ~ rrib, ensures that large logarithms, ln(Mw/mb), are contained in the coefficient 
functions A, B and F only. On the other hand, these functions, which their explicit 
forms can be found elsewhere |7], [8], depend on the top quark mass. For our numerical 
evaluations we use m t = 180 GeV, which is the weighted average of the recent CDF and 
DO results @. 

The terms representing the long distance contributions to the decay rates, enter coef- 
ficients A and B from the one-loop (charm) matrix element of the four-quark operators. 
Besides cc continuum (free quark) contribution, these coefficients receive the pole contri- 
butions from the J/ip and ip', which due to their narrow widths, amount to Breit-Wigner 
terms for these resonances. 

The HME of (1) are calculated by using the trace formalism in the context of HQET 



0: 



< K\v')\sTb\B(v) >= Tr [tf(y')TR(v)M(v,v')\ . (2) 

R l (v') and R(v) are the matrix representations of K % and B respectively, v and v' are 
velocities of initial and final state mesons. M, which represents the light degrees of 
freedom, is related to the Isgur-Wise functions. 

We classify various K- resonances into spin doublets [|T], [|, and thus, using (2), the 
HME to each doublet is expressed in terms of a single universal function. These func- 
tions represent the underlying non-perturbative QCD dynamics, which at present is not 
possible to calculate from first principles and certain model assumptions are required 
for their evaluation. 
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We obtain the following expression for the differential decay rates in the limit of 
massless leptons: 



dT(B -> KH + l~) G 2 F I a 2 



M^y\x + x„) 1/2 \V*V tb \ 2 \Ux)\ 2 



dz 967r 3 y^Tisfy / 

x ((\A\ 2 + \B\ 2 )F l + 2\C\ 2 G l + 2Re[(A + BYC]H i ), (3) 

x = v ■ v , x± = x ±1, y = , z = — 5-, C = s w r, 

where £i(x), (I = C,E,F and G), are the Isgur-Wise functions for each spin doublet 
|], |!|, and q 2 is the invariant mass of the dilepton. The functions F iy Gj and Hi for 
various K-resonances are tabulated in Table 1. For m t = 180 GeV, a s (Mw) = 0.118, 
and 7] = a s (fi)/a s (M w ) for these processes taken to be 1.75, we obtain: 

A = 2.020 + (continuum + resonance terms), 
B = —0.173 + (continuum + resonance terms), 
C = -0.146. 

This shows, excluding the charm loop contributions, the coefficient A is much larger 
than B or C for a large top quark mass. 

Now we turn to the processes B — > K % i/u, which even though are difficult to mea- 
sure experimentally, can be tagged by a large missing energy-momentum. The effective 
Hamiltonian for b — > si/P, arises from the coupling of pure V — A currents, 

Hejf = {^gr) V*V tb [2E(t)](7sL,b)(vL^), (4) 



where t = rn^/mly. The Wilson coefficient E(t) (see ref. ]TT] for the explicit form), is 
not affected by QCD scaling corrections. We take E(t) = 1.66 for m t = 180 GeV. Using 
(2), we obtain the following differential decay rates for B — > K l vv 



dT(B -> K l vv) G\ ( a \ 2 



Mi^^x^^l^^n^^l^E^I^SF,). (5) 



dz 967T 3 \47TS^ / 

The factor 3 mutiplying in (5), is due to the sum over neutrino flavors. 

3 Isgur-Wise function and total decay rates 

In order to estimate the total decay rates, T(B — > K t ££), we need to insert the Isgur-Wise 
functions, in (3) and (5). As we mentioned earlier, in doing so, we need to adopt a 
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model that could also be extended to higher excited final states. The wavefunction model 
of Isgur, Wise, Scora and Grinstein |[2], suits our purpose quite well. The functions £/, 
in this model, can be obtained from the overlap integrals: 



£(W) = \/2L + li L j r 2 dr^* F {r)^ I {r)j l 



Ar 



v.v 



- 1 



(6) 



I and F are initial and final radial wavefunctions respectively, L is the orbital angular 
momentum of the final state meson, and ji is the spherical Bessel function of order L. 
The inertia parameter A, is taken to be fli~3 | 



A 



m s + m q ' 



where m q is the light quark mass. 

We follow IWSG model in using harmonic oscillator wavefunctions with oscillator 
strength f3 for the radial wavefunctions and $p. (3 is assumed to be the same for 
initial and final state mesons in order to satisfy the normalization condition £c(l) — 1 
and £,e,f,g(X) = 0. For example, one obtains 



exp 
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256/3 2 



-m 



K' 



[V.V 



(7) 



for the (0~, 1~) spin doublet by using the ground state radial wavefunction. j3 = 
0. 295GeV is fixed by the best fit of (7) to the experimentally measured B — > (K, K*) + 
(ip, ip') decays ]14|, ||. We use the same value for (3 in our calculation of £f, and^c- 

We would like to remark that eqn. (7) implies the dependence of the universal 
form factor on the K-meson mass. In the heavy quark limit, the members of a spin 
doublet are degenerate in mass, which results in exactly one form factor for each doublet. 
However, for K-mesons, insertion of a realistic value for the mass parameter in (7), leads 
to a different form factor for each member of a doublet. This can be thought of as a 
relaxation of the spin symmetry to some extent which, in our view, results in a more 
reliable estimate of B — > K % transitions. 

The following values are used for our numerical estimates: 



m B = 5.28GeV 
m b = 4.95GeV 



m c - 

\Vts\ 



1.50GeV, 
= 0.043. 



The partial decay widths to various K-resonances are tabulated in Table 2. For B — > 
K l e + e~ and B — > K l fi + fi~, these partial decay rates are different only in the first in- 
terval, ie. small invariant dilepton mass region, where the lower bound of the interval, 
4mf, is determined by the lepton mass. This difference is shown in the first and sec- 
ond columns of Table 2. The enteries of the 3 rd , 4 th , 5 th and 6 th columns of that table 
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re- 



are the partial decay widths for B — > K % £ + £ (£ — e ,//) over ((m^ — 5) 2 , (my, + 5) z 

((m^, + <5) 2 , (m^ - 5) 2 y ({m^ - 5) 2 , (m^ + <5) 2 ), and ((r?v + <5) 2 , (m B - m^*) 2 ) 
spectively. We choose 5 = 0.2 GeV for our numerical estimates. For if* (1680) and 
^(1580), g^ ax = (rriB — rn K t) 2 is in fact smaller than {m^i — 5) 2 and (m^ + 5) 2 respec- 
tively. 

For the decays to the ground state spin doublet (K , K*), we have also examined 
some other parametrizations of the Isgur-Wise function which are commonly used in the 
literature. These are the monopole form, 



u 2 -2 + 2v.v n 

with u = 1.8, and the exponential form, 

£(v.v') = exp [7(1 — v.v')] , 

with 7 = 0.5. The parameter values are obtained from the best fit to the data on 
D -> K£u e . 

In Table 3, we compare the total branching fractions (excluding the resonance con- 
tributions) when different Isgur-Wise functions are used for dileptonic B-decays to K 
and if* (892). The total width of the B is taken to be T(B -> all) = 5.6 x 10~ 13 GeV. 

From (1) and (4), one can calculate the inclusive differential decay rates, ie. 



( a Y 2 \V*V tb 



(B^X£ + D= — — - (1 _ z y 

T(B — > X c eu) dz ^ ) \4ns 2 w ) f(m c /m b ) \V cb \ 2 1 J 



2 



x [(\A\ 2 + \B\ Z )(1 + 2z) + 2|Cf (1 + 2/z) 
+6Re[(A + B)*C}) 



(8) 



and 



dT , „ .. I a • _ 1 ,. 1 ,/, 2 



2 


v t * s v tb \ 


2 


f{m c /m b ) 


v cb \ 


2 



-{B^X a vv) = -— H - . ' * °' |2g(0r x3 



r(B -> X c ez7) cfe v s 7 \^s 2 w ) f(m c /m b ) \V cb \ 2 1 Wl (9) 

x (l - 3^ 2 + 2z 3 ) , 

where 

/(x) = 1 - 8a; 2 + 8x 6 - x 8 - 24a; 4 /n(:r). 

Again, the factor 3 in (9) is due to sum over three neutrino flavors. By normalizing to 
the semileptonic rate in (8) and (9), the strong dependence on the b-quark mass cancels 
out. Using the above expressions, and assuming \V cb \ ~ \V ts \, we obtain: 

T(B -> X s e + e-) = 4.15 x 10~ 18 GeV, 

T(B -> X s n + n~) = 2.81 x 10- 18 GeV, (10) 
F(B -> X a i/P) = 2.43 x 10~ 17 GeV, 
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where the semileptonic branching ratio is taken to be BR(B — > X c eu) = 0.105. 

In Table 4, we have tabulated the total decay rates (excluding the resonance contri- 
butions) and the ratio R = T(B — > K l ££)/T(B — > X s ££) for various K-resonances. 



4 Concluding remarks 

We end this paper with a few concluding remarks. 

From Table 2, we observe that the partial decay widths over the tp resonance region 
(the 3 rd column) are about two orders of magnitude larger than the short distance 
contributions which are mainly from the first region (the 1 st and 2 nd columns for e and \x 
respectively) . The exceptions are the d-wave transitions where this dominance is reduced 
to roughly one order of magnitude. As a result, the total decay rates, T(B — > K l £ + £~), 
are dominated by the long distance (resonance) contributions. The partial decay widths 
over the ip' resonance region (the 5 th column) are not as large, and in fact, except for 
the ground state transitions, the rest are mostly comparable or smaller than the short 
distance contributions. 

The most promising probe of the short distance operators, which are sensitive to 
the new physics, can be done by imposing an upper limit on the invariant mass of the 
dileptons, ie. q 2 (see eqn. (3)), well below m^. As it is reflected in the first and second 
columns of Table 2, for some final state K-resonances, there is a significant difference in 
this partial decay width between e + e~ and pair. This difference is overshadowed in 
the total decay rate by the overwhelming long distance contribution which is independent 
of the dilepton type. On the other hand, much larger partial decay widths is observed 
for some higher excited K-meson states if the above upper cut on q 2 is implemented. For 
example, for q 2 < (m^ - 5) 2 , T(B -> K^(U30)e + e~) and T(B -> #£(1430)// V~) are 
larger than T(B — > K*e + e~) and T(B — > K*[i + fi~) by factors of 6 and 4 respectively 
(see the 1 st and 2 nd columns of Table 2). As we mentioned in the introduction, this is 
not unexpected, as there is a large mass difference between the initial B-meson and the 
final state K-mesons. Taking into account the efficiency for the reconstruction of the 



decaying K-mesons ||15|| , our estimates in this work indicate that along with K*(892) and 
i^(1430), two other P-wave channels, ie. i^i(1270) and i^i(1400), are the best possible 
modes to look for rare dileptonic B-decays in the future B-factories. 

We have compared our results for the ground state K-mesons, K and K*, with those 
obtained from different parametrizations of the Isgur-Wise function. The branching 
fractions, tabulated in Table 3, indicate that, in general, the monopole form of the Isgur- 
Wise function results in a larger estimate of these decay modes. We should also point 
out that the smaller difference between BR(B — > K*e + e~) and BR(B — > K*fi + fi~) 
in IWSG model, is the direct result of the large recoil suppression of the Isgur-Wise 
function in this model. 
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Finally, we note from Table 4, that the exclusive channels we consider in this paper, 
more or less saturate the dileptonic rare B-decay modes. 
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TC^ Nnmp 
ii ± i cl±i±c 


rj 


F l 


rn 




K 


n- 


rp , rp ( 1 _|_ 7/^^ 


,Jb _L «X/ — 




J. \. 1 O JZ< / 


\- 


T i [t* i ( 1 — 7/^^ -I— A. / yy\ 


+ L +^ / V\ ^i'/ ' / J 




J. i. 1 J- i ' J V 7 / 


o+ 




_l «X/ — 


HP 1 T 1 1 — 7/1 


ifi (12701 


i+ 


nr [7* f 1 _|_ <ii\^ -L- A.T 7] 


7> [7* A- A I y>( (*\ — — Z^l 


x [3(l+f/)x + 2Cl — v\xA 


if i (1400) 


i+ 


2/3x_x+ 2 [x+(l - y) 2 +sz] 


2/3x„x+ 2 [x_ + l/z((l - y) 2 x + + z)] 


2/3x_x+ 2 [x_(l + y) +x_] 


K|(1430) 


2+ 


2/3x_x+ 2 [x + (l - y) 2 + 3xz] 


2/3x_x+ 2 [x+ + 3/z((l - yfx + - z)\ 


2/3x_x + 2 [x + (l - y) + 3x_] 


K*(1680) 


l - 


2/3x+x_ 2 [x_(l + y) 2 + xz] 


2/3x + x_ 2 [x+ + l/z{{\ + y) 2 x_ - z)\ 


2/3x + x_ 2 [x+(l - y) + x_] 


K 2 (1580) 


2" 


2/3x+x_ 2 [x_(l + y) 2 + 3xz] 


2/3x + x_ 2 [x_ + 3/z((l + y) 2 x_ + z)] 


2/3x + x_ 2 [x_(l + y) + 3x_] 


if(1460) 


0" 


x+x_(l + y) 2 


x + x_ 


x+x_(l + y) 


K*(1410) 


1" 


x+[x+(l — y) 2 + 4xz] 


x+[x++4/z(((l+y) 2 x_ + z)] 


x+[3(l-y)x+ + 2(l + y)x_] 



Table 1: The functions Fj, Gj and i7j for various K-resonances. 



Name 


4m 2 , (m^ — 5) 2 


4m 2 , (rrty - <5) 2 


(m^, - <5) 2 , (rriip + S) 2 


(m^p + 5) 2 , (m^> - 8) 2 


(m^ -8) 2 ,( mi ,, +8) 2 






1.1 x l(T iy 


1.1 x 10" i9 


3.3 x 10- iY 


6.2 X 10" 2U 


2.4 x 10~ i8 


5.6 x 10^ 2U 


iT* (892) 


2.1 x l(T i9 


1.9 x 10" i9 


8.0 x 10- iY 


1.8 X 10~ i9 


8.3 x 10" i8 


1.8 x 10~ iy 


K*(1430) 


8.1 x 10~ 2U 


8.0 x 10- 2U 


9.1 x 10- iS 


6.6 x 10- 2i 


2.7 x 10~ 2U 


3.2 x 10~ 23 


#i(1270) 


2.9 x 1(T 19 


1.8 x 10~ 19 


2.7 x lO" 17 


2.4 x 10~ 2U 


2.6 x 10~ iy 


7.3 x 10" 22 


if 1 (1400) 


7.2 x 10~ iy 


5.7 x 10" 19 


5.0 x 10" 1Y 


3.8 x 10~ 2U 


1.6 x 10" iy 


2.3 x 10~ 22 


#2(1430) 


1.2 x 1(T 18 


7.5 x 10~ 19 


7.2 x 10- 17 


5.3 x 10" 2U 


2.0 x 10~ iy 


2.3 x 10~ 22 


K*(1680) 


3.3 x 1(T 2U 


2.1 x 10" 2U 


1.7 x 10~ 19 


5.8 x 10~ 23 






#2(1580) 


1.1 x 10" 19 


5.7 x 10~ 2U 


1.1 x 10" 18 


4.4 x 10" 22 







#(1460) 


6.7 x 1(T 2U 


6.7 x 10~ 2U 


5.0 x 10~ iS 


2.8 x 10~ 2i 


2.2 x 10~ 2i 


1.3 x 10~ 24 


IT* (1410) 


2.5 x 10" 19 


1.8 x 10" 19 


3.1 x 10" 1Y 


2.5 x 10~ 2U 


1.3 x 10 -19 


1.7 x 10" 22 



Table 2: Partial decay widths (in GeV) for B — > K l £ + £ (£ = e, jj). The difference between electron and muon 



final states is in the first interval only (small invariant mass for dilepton pair) which is shown in 
the first two columns. 
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Isgur — Wise 
function 


B -► i-Ce+e" 


B -► 


B -► Kz/P 


B -► K*e+e" 


B -► K*fi+fi- 




Harmonic 
wavefunction 
(3 = 0.295 


3.4 x 10~ 7 


3.4 x 10~ 7 


2.9 x 10~ 6 


8.9 x 10~ 7 


8.6 x 10~ 7 


8.0 x 10" 6 


Monopole 
u = 1.8 


6.6 x 10~ 7 


6.6 x lO" 7 


5.7 x 10~ 6 


3.8 x lO" 6 


2.5 x 10~ 6 


2.1 x 10~ 5 


Exponential 
7 = 0.5 


2.9 x 10~ 7 


2.9 x 10~ 7 


2.5 x 10" 6 


3.0 x 10" 6 


2.1 x 10~ 6 


2.0 x 10~ 5 



Table 3: Comparison between the branching fractions of the rare B-decays to the ground state K-mesons 
(resonance contributions are excluded) for different Isgur- Wise functions. 



K l Name 




Mass (MeV) 


B -► i^ l e-e+ 


B -► K l fi-j2 + 


B -► K l vv 


K 


o - 


497.67 ± 0.03 


1.9 x 10~ iy , (4.6) 


1.9 x 10~ iy , (6.8) 


1.6 x 10~ 18 , (6.6) 


K*(892) 


1- 


896.1 ±0.3 


5.0 x 10" iy , (12.0) 


4.8 x 10" iy , (17.1) 


4.5 x 10~ 18 , (18.5) 


K*(1430) 


0+ 


1429 ± 7 


6.5 x 10~ 2U , (1.4) 


6.5 x 10~ 2U , (2.3) 


6.0 x 10~ iy , (2.5) 


#i(1270) 


1+ 


1270 ± 10 


2.7 x 10~ iy , (6.5) 


1.6 x lO -19 , (5.7) 


1.3 x 10~ i8 , (5.3) 


#i(1400) 


1+ 


1402 ± 7 


6.3 x lO" 19 , (15.2) 


4.8 x 10~ iy , (17.1) 


3.6 x 10~ i8 , (14.8) 


K|(1430) 


2+ 


1425.4 ± 1.3 


1.1 x 10~ 18 , (26.5) 


6.3 x 10~ 19 , (22.4) 


4.5 x 10~ 18 , (18.5) 


K*(1680) 


1- 


1714 ± 20 


3.1 x 10~ 2U , (0.7) 


1.9 x 10~ 2U , (0.7) 


1.5 x 10~ iy , (0.6) 


K 2 (1580) 


2 - 


» 1580 


1.0 x 10~ iy , (2.4) 


5.1 x 10~ 2U , (1.8) 


3.0 x 10" iy , (1.2) 


K(1460) 


o- 


« 1460 


5.4 x 10~ 2U , (1.3) 


5.4 x 10~ 2U , (1.9) 


4.7 x lO" 19 , (1.9) 


K*(U10) 


1- 


1412 ± 12 


2.2 x 10~ iy , (5.3) 


1.5 x 10~ iy , (5.3) 


1.3 x 10~ 18 , (5.3) 




Sum: 76.1% 


Sum: 81.1% 


Sum: 75.2% 



Table 4: The total decay rates, in GeV (resonance contributions excluded), and the ratio 
R = F(B -> KH1)/T{B -► X s ll) (in %) for various K-resonances. 
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